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ABSTRACT 



We present 32 epochs of optical (3300-9700 A) spectrophotometric observations of the nearby quintessential "normal" type la super- 
nova (SN la) SN 201 lfe in the galaxy M101, extending from - 15 to + 97 d with respect to S-band maximum, obtained by the Nearby 
Supernova Factory collaboration. SN 201 lfe is the closest (ji = 29.04) and brightest (S max = 9.94 mag) SN la observed since the 
advent of modern large scale programs for the intensive periodic followup of supemovae. Both synthetic light curve measurements 
and spectral feature analysis attest to the normality of SN 201 lfe. There is very little evidence for reddening in its host galaxy. The 
homogeneous calibration, intensive time sampling, and high signal-to-noise ratio of the data set make it unique. Thus it is ideal for 
studying the physics of SN la explosions in detail, and for furthering the use of SNe la as standardizable candles for cosmology. Sev- 
eral such applications are shown, from the creation of a bolometric light curve and measurement of the 56 Ni mass, to the simulation 
i—i of detection thresholds for unburned carbon, direct comparisons with other SNe la, and existing spectral templates. 
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(t$ 1. Introduction uvm2, and uvwl filters. BVRI photometry obtained over six 



months appeared in Richmond & Smith ( 20 12}, Vinko et al. 



Having exploded in the Pinwheel Galaxy just 6.4 Mpc distant ( |2012| l and |Munari et al.| ( |2013| l. |Tammann & Reindl| ( |201 1\ 
(z = 0.00080 + 0.00001 ; fi = 29.04 ± 0. 19, |Paturel et al.]2003l used extensive optical photometry of SN 201 lfe obtained by the 
Shap pee & Stanek|2011] ), the Type la supernova SN 201 lfe rep- American Association of Variable Star Observers (AAVSO) and 

the tip of the red-giant branch distance to its host (M101) to mea- 
sure Hq. Covering the near-infrared, Matheson et al. (2012) pre- 



resents a rare opportunity for intensive study. Its discovery (Nu 



gent et al.|2011a[) by the Palomar Transient Factory (PTF; 



.aw 



et al.||2009[ |Rau et al.|[2009) less than 12 hours after outburst sented high-cadence JHK photometry from the Wisconsin Indi- 



( Nugent et al.||2011bi precipitated a number of diverse ground- ana Yale NOAO (WIYN) telescope with WHIRC. The apparent 



based and space-based follow-up campaigns. Promptly initiated, brightness (V ~ 10 at peak) of SN 201 lfe made multi-epoch 
high-cadence observations from a number of these have been spectropolarimetry much more accessible than usual ([Smith 



published already. |Brown et al.| ( 2012] > presented two months |et al.|[20TT) . Finally, |Parrent et al.| ( |2012| have published 18 



of nearly continuous Swift/UVOT follow-up in ultraviolet uvw2, optical spectra from multiple telescopes, starting 1.2 days after 
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explosion with a 1.8 day average cadence, along with a spectro- 
scopic analysis including constraints on unburned carbon. 

These data sets and others can address numerous long- 
standing questions about the nature of the Type la supernova 
(SN la) progenitor systems, environs, and explosion mecha- 
nisms. It has generally been thought that SNe la arise from the 
thermonuclear disruption of white dwarf stars accreting material 
from a companion donor star (for an overview, see e.g. Branch 
T995J. Using images obtained just hours after outburst, 



et al. 



pre-explosion X-ray limits, and the inferred Ni yield, 



Bloom 



et al. (2012 1 constrain the progenitor primary of SN 201 lfe to 
be a white dwarf or neutron star - the first direct imaging ev- 
idence for a compact primary. Pre-explosion multi-wavelength 
archival images of the SN 201 lfe stellar neighborhood require 
the mass-donating secondary to be another white dwarf, sub- 
giant, or main-sequence star; red giants and helium stars are ex- 
cluded with the companion restricted to M < 3.5 M G ( [Li et al.| 
201 1)>. Radio and x-ray observations probing the circumstellar 



environment of SN 201 lfe imply a progenitor system with a 
mass loss rate as low as 6 x 10~ 10 M yr~', severely constraining 
models where the secondary is not a white dwarf ( |Chomiuk et al.| 
|2012[|Horesh et al.|20"l2[|Margutfi et al.|2012) >. A case study in 
attempting to constrain explosion models from spectrophotomet- 
ric observations slightly favors a white dwarf companion (Rtipke 



|et al.|2012| > but further detailed modeling is needed to be more 
conclusive. The future emergence (or absence) of a surviving 
non-compact secondary, brightened as a consequence of its in- 
teraction with the SN ejecta, may be decisive in identifying the 
progenitor system (Shapp ee et al.|2012] >. 

Such progress is exciting for several reasons; chief among 
these is the potential ability for such results to reduce uncertainty 
about the reliability of SNe la as tools for observational cos- 
mology. The use of SNe la as standardizable candles (Phillips 
|1993] > brought about the discovery of the accelerating expansion 
of the Universe just over a decade ago ( jRiess et al.|1998[ [Perl 



mutte r et al.|1999| l. SNe la have since become a key means for 



constraining the unknown physics of cosmic acceleration, called 
"Dark Energy" (e.g., |Guy et aL|2010[|Howell|201 1 [[Suzuki et al. 
2012). A fundamental physics result is thus inextricably linked 
to the details of binary stellar evolution and stellar death, details 
which SN 201 lfe may help sort out. It is thus fortuitous that 
SN 201 lfe is also a normal SN la, and is highly representative 
of the typical SN la sought for placement on a Hubble diagram. 

In this article, the Nearby Supernova Factory (SNfactory, 
Aldering et al. 2002) presents an atlas of 32 spectrophotomet- 
ric observations of SN 201 lfe extending from - 15 to +97 d 
with respect to the time of maximum light. Some of these spec- 
tra have appeared in the Ropke et al. (2012) study and a quick- 
pipeline reduction of the first spectrum presented here has ap- 
( 2012[ ). This atlas should become a use- 



peared in Parrent et al. 



ful resource in studying SN la physics, in exploring systematics 
in the analysis of SN la spectra and spectral indicators, and in 
constructing spectral templates for SN cosmology applications. 

The remainder of this article is organized as follows. In 
Sect. [2] we describe the observations and data reduction proce- 
dure, present the reduced spectrophotometric time series, syn- 
thesize photometry from it, and analyze the synthetic light 
curves using a standard SN la light curve fitter. Analyses of 
light curve residuals, possible extinction in the host galaxy, and 
the bolometric light curve appear in Sect. [3] This section also 
covers spectral feature measurements, sub-classification within 
existing schemes, parameterized spectral fitting, and examina- 
tion of unburned carbon signatures. In Sect. [4] we place our 
data and analysis in context but primarily seek to demonstrate 



some useful features and applications of the data set. Conclu- 
sions appear in Sect. [5] The spectrophotometry of SN 201 lfe 
is available for download in electronic form at the SNfactory 
project websiteQand at the CDS via anonymous ftp to cdsarc.u- 
strasbg.fr (130.79.128.5) or via |http : //cdsweb . u- str a sbg . 
|fr/cgi-bin/qcat?J/A+A/| 



2. Spectrophotometric observations 

2.1. Data acquisition 

The data were obtained using the SuperNova Integral Field Spec- 
trograph (snifs, |Lantz et al.|2004| ). snifs is a fully integrated in- 
strument optimized for automated observation of point sources 
on a structured background over the full ground-based optical 
window at moderate spectral resolution. It consists of a high- 
throughput wide-band pure-lenslet integral field spectrograph 
(IFS, "a la TIGER;" [Bacon et al.|1995l[2U0T) , a multi-filter pho- 
tometric channel to image the stars in the vicinity of the IFS 
field-of-view (FOV) to monitor atmospheric transmission dur- 
ing spectroscopic exposures, and an acquisition/guiding chan- 
nel. The IFS possesses a fully-filled 6'.'4 x 6'.'4 spectroscopic 
field of view subdivided into a grid of 15 x 15 spatial elements, 
a dual-channel spectrograph covering 3200-5200 A and 5100- 
10000 A simultaneously with FWHM resolutions of 5.65 A and 
7.54 A respectively, and an internal calibration unit (continuum 
and arc lamps), snifs is continuously mounted on the south 
bent Cassegrain port of the University of Hawaii 2.2 m telescope 
(UH 2.2m) on Mauna Kea. The telescope and instrument, under 
script control, are supervised remotely. 

PTF discovered SN 201 lfe (PTFllkly) in images obtained 
on 201 1 Aug. 24.2 (UTC, used throughout). |Nugent et aL 
(2011b|) derive an explosion date of Aug. 23.7. SNfactory 



follow-up observations commenced on Aug. 26.3, an estimated 
2.6 days after the explosion, and encompass 32 nights with two 
consecutive spectrophotometric exposures on most of the nights. 
The observing log is shown in Table [T] Daily observational ca- 
dence was maintained until 12.5 days after explosion. Lack of 
telescope access then forced a four day gap, but the daily ca- 
dence was re-established for days 16.5 to 21.5 after explosion. 
The cadence was relaxed to every 2-3 days until SN 201 lfe was 
no longer observable from Hawaii, 41.5 days after explosion. 
Follow-up resumed a month and a half later, when observations 
were made over 7 nights spanning 20 days. The last spectra re- 
ported here were obtained Dec. 16.6, about 115 days after explo- 
sion. The position of SN 201 lfe on the sky during the follow-up 
campaign forced observations at low altitude. The mean airmass 
was 2.1 and most of the early observations were taken during 
astronomical twilight. 

2.2. Data calibration 

All spectra were reduced using SNfactory's dedicated data re- 
duction pipeline, similar to that presented in Sect. 4 of |Bacon] 
|et al.| ( |2001[ ). A brief discussion of the software pipeline is 
presented in Aldering et al.| ( |2006[ ) and is updated in Sect. 2.1 
of |Scalzo et al.| ( |2010| ). SN201 lfe is an unusual target for snifs 
and its associated data reduction software pipeline. The bright- 
ness of the SN itself and the small amount of time that the target 
was above the horizon imposed constraints on integration time. 
This affects our ability to use our standard multi-filter photo- 
metric channel for flux calibration on non-photometric nights. 



http : // snf actory . lbl . gov 
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R. Pereira & the SNfactory: Spectrophotometric time series of SN 201 lfe 
Table 1. Observing log for snifs spectra of SN 201 lfe 
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Notes. In the case of multiple consecutive observations in a single night, the phases, dates and airmass correspond to the middle of the first 
exposure, while the seeing is the average value. 

(a) Days relative to the date of explosion derived by Nugen t et al. ' 2011b i: MJD 55796.696. 

(b) Phase, observer-frame days relative to S-band maximum light: MJD 55814.51. 

(c) JD - 2400000.5 

(d) • photometric, o non-photometric, (#) number of standard stars observed during the night and used for atmospheric extinction and telluric 
absorption correction. 



Also, some observations were obtained using interrupt time lim- 
ited to ~ 1 h which did not allow for execution of more than 
one standard star observation needed for optimal calibration. To- 
gether, these factors limit the accuracy of our flux calibration of 
SN 201 lfe, despite the brightness of the target. We will thus 
elaborate on the data reduction discussion given in |Scalzo et al.| 
(2010) and on modifications needed in order to flux calibrate the 
observations presented here. 

Following standard low-level processing, CCD frames from 
both snifs channels are mapped into (x, y, A) datacubes. Spec- 
tra are extracted from these using a chromatic spatial point- 
spread function (PSF) assuming spatially flat background, thus 
removing most of the host galaxy light as part of the sky. Host 
galaxy light contamination was evaluated using a background- 
subtracted stack of 5 x 60 s CFH12k /?-band exposures of the 
host galaxy of SN 201 lfe (M101), which was photometrically 
calibrated using SDSS r-band observations of the same region. 
Several snifs acquisition images were stacked and astrometri- 
cally aligned to the CFH12k stack in order to determine the lo- 
cation of SN 201 lfe within the CFHT image. The mean r-band 
host galaxy surface brightness, inside the snifs FOV centered 
around the SN, equals 22.15 mag arcsec~ 2 and is almost flat, 
with structured residuals in the order of r ~ 23 mag inside 1" 



diameter apertures. The remaining background structural effects 
are thus considered negligible with respect to the brightness of 
SN 201 lfe, so the detailed host galaxy subtraction method de- 
scribed by Bongard et al. ( 201 1) was not necessary. 

The extracted spectra are next merged and truncated to the 
final wavelength range of 3300-9700 A. Flux calibration is per- 
formed using an instrumental flux solution derived from all spec- 
trophotometric standard stars observed during the same night, 
and a mean atmospheric extinction computed from 6 years of 
snifs data (But on et al.|2013[ ). Observations of multiple standard 
stars were not always possible, particularly on interrupt nights. 
This limits a proper nightly extinction computation, so for the 
sake of consistency we use the mean extinction law on all nights. 
In this particular framework, the precision of the flux calibration 
depends on the number of standard stars observed per night, but 
not on their airmass distribution. 

During non-photometric nights, any achromatic (grey) dif- 
ferential atmospheric attenuation between the observations of 
SN 201 lfe and the standard stars is accounted for using simul- 
taneous observations of field stars, through the snifs multi-filter 
photometric channel. Because SN 201 lfe was very bright, rela- 
tively short exposure times were used to avoid saturation of the 
spectrograph, so the number of field stars visible at all epochs 
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under very different attenuation conditions is small. After de- 
tection of the objects in the multi-filter field, the object cata- 
logs were manually inspected to insure proper star selection and 
to prevent inaccurate inter-epoch astrometric and photometric 
alignments, which are performed using a modified version of 
the SuperNova Legacy Survey photometry code (poloka, Astier 
et al.|2006) . The flux ratio between each exposure and a refer- 
ence exposure (with the best seeing) equals the integral of the 
convolution kernel needed to degrade the reference image PSF 
into the "photometric frame" of the other image. This night-to- 
night variation in the brightness of the field stars, normalized 
by observations on photometric nights, provides the photomet- 
ric scaling between fields and is used as a per-exposure grey 
flux correction to the mean atmospheric extinction, thus allow- 
ing absolute flux calibration. Seeing and night photometricity 
were assessed using quantitative analyses of snifs guider video 
frames acquired during our exposures, along with deglitched 
CFHT Skyprobe data (Cuillandre et al. 2002 1, attenuation es- 
timates from the instrumental flux solution and photometric ra- 
tio computations (cf. Sect. 5 of B uton et al.|2013[ ), and infrared 
satellite imagery. 

The flux-calibrated spectra were corrected for telluric ab- 
sorption using a nightly average spectrum determined using 
standard star observations. While consistent with the errors, 
this correction is nevertheless imperfect, and some (very small) 
glitches are visible in the spectra, especially for later phases at 
A > 9000 A. Milky Way dust extinction along the line of sight 
is corrected assuming E(B - V) — 0.0088 mag (Schlegel et al. 
1998|l and an exti nction law with R v = 3.1 ( |Cardelli et al.|1989| 
O'Donnell 1994). No extinction correction for dust in the host 



galaxy is applied (see Sect. 3.2 1. Multiple spectra obtained the 



same night were merged using a variance-weighted mean. 



2.3. Spectral time series and synthesized photometry 

The final spectrophotometric time series of SN 201 lfe is de- 
picted in Fig. [T] The time series possesses a number of out- 
standing qualities, beyond the extremely high S/N, despite the 
calibration challenges. Each spectrum spans the entire ground- 
based optical wavelength window. Excluding an inconvenient 
four day gap, daily sampling is achieved on the rising side of 
the light curve. Even though SN 201 lfe observations were per- 
formed on non-standard conditions when compared to the typi- 
cal snifs target, our observing strategy and calibration procedure 
do an excellent job of removing atmospheric artifacts. These 
features, combined with the cosmological utility of SN 201 lfe 
as a "normal" SN la, should make these data a highly useful re- 
source. 

Synthetic light curves were generated by integrating the 
product of each spectrum with a set of non-overlapping top- 
hat filters with perfect photon transmission in the following 
wavelength ranges: 3300-4102 A, 4102-5100 A, 5200-6289 A, 
6289-7607 A, and 7607-9200 A, respectively UBVRI SN f This 
filter set avoids the split between both snifs spectrograph chan- 
nels but is contiguous otherwise. For convenience these mea- 
surements are given in the Vega magnitude system by apply- 
ing zero-points computed from the latest Hubble Space Tele- 
scope (HST) spectral observation of Vega (Boh lin|2007[ ), avail- 
able from the CALSPEC database]^] The synthesized photometry 
is listed in Table[2j along with the integrated flux over the whole 
snifs optical window. For ease of comparison with other SNe, 



we also present synthetic light curves created using the latest 
photonic responses and zero-points given by Bessell & Murphy 
( |20T2l BM12), which are denoted as 5V/?/ B Mi2 o in Table|2j After 
interpolation of those filter responses into 10 A steps using cu- 
bic splines, our integrator agrees at the 0.001 mag level with the 
colors measured by Bessell & Murphy (2012) using the afore- 
mentioned Vega spectrum and the zero-points derived by these 
authors. 

The achromatic absolute flux calibration accuracy for each 
night as estimated by the calibration pipeline is shown in the 
cr column of Table |2] This value takes into account the errors of 
the spectral extraction using the chromatic spatial PSF, estimated 
from the mean scatter of flux calibrated standard star observa- 
tions during all the photometric nights of the SNfactory data set, 
as well as the errors of the instrumental flux solution, the mean 
atmospheric extinction and the photometric ratio computations. 
The main contributions are the achromatic extraction uncertainty 
of 0.03 mag and 0.02 mag empirically found for bright and faint 
standard stars (Buton 2009), and a 0.03 mag additional uncer- 
tainty on non-photometric nights due to the photometric ratios 
usage ( |Pereira||2008[ l. The chromatic spatial PSF flux extrac- 
tion efficiency was evaluated on residual spectra, extracted from 
the PSF-subtracted datacubes by performing aperture photom- 
etry centered at the position of the SN and within a radius of 
3cr seeing. The median missing flux for all exposures, on the 
5 broad-band filters of the SNfactory filter set, was found to be 
below the half-percent level for all bands except for t/sNf > where 
it is ~ 0.7%. A & 1% residual color effect is also discernible, 
when comparing Z?sni with the other bands at high airmasses 
(> 2.7), and is due to the fact that at such atypically large air- 
mass, the atmospheric differential refraction (ADR) is so large 
that too much light is projected out of the snifs FOV The con- 
tribution to the total error from both residual chromatic effects is 
much smaller than that from the absolute flux calibration. 



2.4. SALT2 fit 

The UBVRsm light curves were simultaneously fitted using 
SALT2.2 ( |Guy et al-pOTO) SALT2 henceforth), for estimation 
of the date of maximum light and assigning phases to the ob- 
servations (relative to maximum light). The publicly available 
SALT2 spectral mode^jseems to attain the (synthesized) B peak 
slightly before phase 0, which is a nuisance when we use it for a 
precise determination of the date of maximum light of our light 
curve. The shift needed to be applied to the fitted date of maxi- 
mum, in order for it to match with phase in B, was determined 
empirically to be ADayMax B * -0.7 d. The SALT2 code used 
here was modified to take this shift into account. 

For a coherent treatment of the SNfactory filter set by 
SALT2, which uses a BD+17 4708 based magnitude system, the 
following zero-points are added to the BD17-snls3.dat mag- 
nitude system configuration file: 9.787, 9.791, 9.353, 9.011, 
and 8.768; which represent BD+17 4708's synthetic magnitudes 
with relation to the Vega spectrum, in UBVRIsNi respectively. 
Equivalently, if we intend to use SALT2 to fit light curves syn- 
thesized using the UBVRIbmii photonic responses, the zero- 
points to be added to the configuration file are 9.709, 9.902, 
9.469, 9.163, and 8.843. It should be stressed that every SALT2 
fit result presented here is in the "standard" SALT2 UBVRI mag- 
nitude system, allowing direct comparison with the literature. 
This magnitude system is based on the Bessell ( 1990 ) filter trans- 



2 |ftp:// ftp . stsci . edu/cdbs/current_calspec/alpha _lyr_ 
stis_8Q5.fits 



] 



3 http : // supernovae . in2p3 . f r/~guy/salt/download_ 
templates .html 
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Fig. 1. snifs spectrophotometric time series of SN 201 lfe from -15 to 100 days relative to S-band maximum light. Breaks in the axis on the 
right indicate gaps and changes to the observing cadence. The first break corresponds to a four-day gap in daily cadence before maximum. The 
second marks the change from daily to alternating two/three day cadence. The final break is a 50 day hiatus imposed by lack of accessibility to the 
target from Mauna Kea. 
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Fig. 2. SN 201 lfe synthesized light curves using the UBVRIsm filter set. Filled and open symbols stand for photometric and non-photometric 
nights respectively. The results of a SALT2 simultaneous fit of UBVRsw in the phase range - 16 < / < + 25 d are shown as solid lines, along with 
the corresponding residuals (SALT2 - snifs) on the lower panel. The shaded areas represent the SALT2 model error. The residuals for the first 
points of UsNf and B SNf fall outside the panel, and the RMS on the residuals for each band ignores the first 2 points. The break in the time axis 
corresponds to the ~ 50 day gap in follow-up during which SN 201 lfe was not visible during the night from Hawaii. Note the change of scale of 
the extended time axis covering the late observations. 



missions, shifted to match Landok] (JT992J observations. If we 
use instead the magnitude system based on the BM12 transmis- 
sions, the effect on the fitted parameters is very small, and com- 
pletely within the quoted uncertainties. 

The fit was performed using all available observations within 
four weeks of maximum light, including the very first one, which 
falls just outside the default SALT2 phase range of - 15 < t < 
+ 45 d. Early observations are deweighted by the S ALT2 error 
model. Simply omitting the first four observations (phases < 
- 12 d) results in a shift in the date of maximum of -0.018 d, 
well inside the quoted error. 

The light curves, fits and corresponding residuals can be 
seen in Fig. [2] (for a comparison with light curve results already 
published on the literature see Sect. |3.1| l. The fitted B maxi- 
mum of 9.94 ± 0.01 mag was reached on 201 1 September 10.5. 
The best-fit SALT2 parameters are x x = -0.206 + 0.071, c = 
-0.066 + 0.021 (see Table [3] for a full summary of extracted 
photometric and spectroscopic parameters). Excluding f/SNf- 
band data from the fit, we obtain a negligible shift for the date 



and magnitude at maximum, and X\ = -0.149 ± 0.096, c = 
-0.061 ± 0.027. The light curve shape parameters are typical 
of a "normal" (if slightly blue) SN la: the median values for the 
x\ and c distributions of the nearby (z < 0.1) SNe la used by 



Conley et al^pOTl) are respectively -0.249 and -0.026 (J. Guy, 
private communication). The V-band absolute magnitude (as- 
suming fi = 29.04 + 0.19, |Shappee & Stanek]|2011Q at B peak, 



S V = -19.05 + 0.19 mag matches perfectly the average 



found by|Riess et al] ([20091 -19.06 + 0.05 mag). 



The RMS of the residuals of the four fitted filters, ignoring 
the first two exposures due to the large discrepancies with the 



SALT2 model (cf. Sect. |4\2) , are respectively 0.06, 0.03, 0.04 
and 0.04 mag. The points that deviate the most from the fit 
are those for phases t < -10 d, showing the inadequacy of the 
SALT2 model for such early phases. Nevertheless, the exposure 

of night MJD 55805 (t 9 d, cf. Fig. |2| seems to present a 

systematic error in its flux calibration of ~ 0.1 mag, since it is 
brighter than we would expect based on observations on neigh- 
boring nights. 
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3. Analysis 

3. 1 . Light curve comparison 

The optical light curves of SN 201 lfe published by different 
authors allow us to assess the reciprocal precision of multiple 
independent observations of this SN, and ultimately the bene- 
fits of using SNe la spectrophotometry for synthetic light curve 
creation. We use for these comparisons the first set of SNfac- 
tory observations, up to MJD 55840. For every comparison light 
curve provided by other followup campaigns, the approach used 
to compute the residuals with respect to the reference snifs is 
the same: a light curve is synthesized from our spectrophoto- 
metric time series in the same passband as the comparison light 
curve, and then fit using gaussian processes (as implemented 
by scikit-learn; Pedregos a et al.|2011| l employing a squared 
euclidean correlation model weighted by measurement errors. 
The fitted model is then used to evaluate the synthesized snifs 
light curve at the same observation dates of the comparison light 
curves. The magnitude residuals are always reference (snifs) mi- 
nus comparison. Positive or negative residuals thus mean that 
the comparison observation is respectively brighter or fainter 
than the reference. 

We start with a brief comparison with the space-based 
Swift/UVOT observations by |Brown efaL] pOTI} B12), before 
performing an extensive comparison with the three ground-based 
photometric followup compaigns who published observations 
of SN 201 lfe: |Richmond & Smith| (|20T2l RS 12), |Vink6 et al.| 
pOTIl V12) and |Munari et al.| ( |2013| Ml 3). This trio of light 
curves and residuals with relation to synthesized snifs photome- 
try using BM12 passbands and zero points, are shown in the two 
upper panels of Fig. [3] While a quick look at the superposed light 
curves may give the impression of a good agreement between all 
experiments, a careful study of the residuals shows that this is 
not always the case. The median and normalized median abso- 
lute deviation (nMAD) of the residuals between all the experi- 
ments are summarized in the matrices shown in the lower two 
panels of Fig. [3] These matrices relate the statistics of the resid- 
uals (observational bias and scatter) between a specific data set 
and all the other ones, where each row uses a common reference 
and each column represents a different comparison data set. The 
bias and scatter matrices are not perfectly (anti-)symmetric since 
that the interpolation of the light curve depends on the time sam- 
pling and measurement errors of the reference data set, allied to 
the fact that the nMAD is less dependent on outliers. The statis- 
tics of the residuals with relation to snifs, plotted on the second 
panel of Fig. [3] can be read from the topmost row of each matrix. 



3.1 .1 . |Brown~eTan i |20T2) 



B12 present measurements in several UVOT filters, of which 
only b and v ( |Poole et al.| [2008 ) fully overlap with the snifs 
wavelength range. Unfortunately both filters saturated before 
maximum light, and as such our comparison is limited to phases 
t < - 8 d and + 18 < t < + 25 d. We use the most recent UV OT 
effective area curves and zero points ( |Breeveld et al.|20TT) and 
integrate our spectra in F v . This approach was validated us- 
ing a Vega spectrum, from which we found exactly the same 
zero points as reported by |Breeveld et al.| ( |2011| l. We found 
a clear bias on the residuals, whose median values are -0.14 
and -0.20 mag respectively for b (with comparison dates only 
at t < - 8 d) and v, with no evident correlation with magni- 
tude. This bias also seems to be present for the bluer filters (cf. 
Sect.|3.3[). The cause for this effect is not clear. The photometric 



calibration of UVOT is a notoriously difficult endeavor (Poole 
|et al.|2008] l, and the fact that the instrument nears the saturation 
limit for SN 201 lfe could point to problems in the coincidence 
loss or aperture corrections. S -corrections ( |Suntzeff[ 2000) may 
also play a role, since |Brown et al. ( 20 1 2) > use the |Hsiao et al. 



( 2007) > template to convert the observed count rates into flux, 



while these template spectra do not exactly match the observed 
ones at early phases (cf. Sect.|4.2[). 



3.1 .2. [Richmond & Smith] |2012) 

RS12 present BVRI observations coming primarily from the 
Rochester Institute of Technology (RIT) observatory. These 
observations are the most extensive and consistently calibrated 
from the full AAVSO followup sample used by Tammann & 
|Reindl| ( |201 1) . The biases for VRI are small and within the 
quoted errors of the synthetic photometry presented in this work, 
while the scatters are smaller or similar to the ones found for the 
SALT2 fit. The B-band however shows a large bias with rela- 
tion to snifs observations (~ 0.1 mag), due to an S -correction 
problem with the calibration of RS12. That is evidenced by the 
unusually large color coefficient derived by the authors for this 
filter (0.24 ± 0.04), and by the fact that when using the RS12 
B-band light curve as reference (second row of the first matrix), 
we find large biases with relation to every other data set. The 
CCD used in the RIT observatory has a very pronounced quan- 
tum efficiency (QE) drop in the B-band wavelength range, which 
explains this effect. Using CALSPEC spectra for stars with sim- 
ilar colors to the ones from the PG 163 3 +099 standard field used 
by RS12 for their calibration, the BM12 passbands convolved 
by the KAF-1600 CCD QE curve, and an atmospheric extinc- 
tion model (pyExtinction; But on et al.|2013 1 based on typical 
atmospheric values for the RIT site, we obtain a b - v color co- 
efficient of 0.31 ±0.10, while that ignoring the QE effect gives 
-0.07 ±0.08. The QE convolved passbands also improve the cor- 
relation of the B-band synthetic photometry with the instrumen- 
tal Z?sn _ ^Vega> after inversion of the color correction performed 
by RS12. Differential ADR between a blue SN la and red ref- 
erence stars will also affect filter photometry at some level, es- 
pecially at the high airmasses of SN 201 lfe observations. The 
authors' derived values for the B-band date of maximum, peak 
magnitude and decline rate are then naturally different from ours, 
being greater by approximately 1 day, 0.05 mag and OTFlrespec- 
tively. We point out that the same values derived by |Tammann & 



|Reindl| ( |201 1 1 using the full AAVSO sample are perfecly com- 
patible with those found in the present work. 

3J.3. |Vink6eraLl(20T2) 

V12 published BVRI photometry obtained at the Konkoly Ob- 
servatory. The BI bands are more compatible with snifs obser- 
vations than their RS12 counterparts, but this is not the case for 
VR. The CCD that was used has a less steep QE curve in the bluer 
wavelengths, and is thus less prone to S -correction problems as 
shown by the smaller bias for the B-band, still ~ 0.04 mag fainter 
than the snifs photometry. This bias does not seem to be present 
if we use M13 (cf. |3.1.4[ ) instead as the reference light curve, 
which leads us to think that the bias is affected by the quality 
of the interpolated light curve, and its dependence on individ- 
ual point measurement errors. Globally, V12 data shows some 



4 IVinko et al.| 42012b attribute this difference to a probable misprint, 
and find Am ls (B) = 1.12 ± 0.05 using RS12 data. Fitting the same data 
with SALT2 we find 1.17 ± 0.04. 



Article number, page 8 of 



23 



R. Pereira & the SNfactory: Spectrophotometric time series of SN 201 lfe 




Fig. 3. First panel: comparison of synthesized BVRI snifs light curves of SN 201 lfe, using BM12 passbands and zero points, with published 
photometry from |Richmond & Smith] J2012| RS12), |Vm"k6 et al.|pOT2l V12) and |Munari et al.lpOB) M13). Second panel: residuals between 
all light curves and the snifs interpolated light curve, whose uncertainty is represented as the gray shaded region. Third & forth panel: median 
and normalized median absolute deviation of the residuals for each combination of reference (row) and comparison (column) light curves. The 
statistics of the residuals shown in the second panel, with snifs as reference, are on the topmost row of each matrix. The RqIq filters used by Ml 3 
are not directly comparable to RI from RS12 and V12. 



tension when used as reference, especially with snifs for BVR. 
The scatter of the residuals with relation to snifs is of the order 
of 0.05 mag for all bands, which is higher than RSI 2, especially 
for the redder filters. The Z?-band date of maximum and decline 
rate obtained by V12 using MLCS2K2 are in accord with our 
own. 



3.1 .4. [Munari et al.|j2013) 



Finally, Ml 3 published BVRqIq photometry obtained using sev- 
eral telescopes of the Asiago Novae and Symbiotic stars collab- 
oration (ANS; Munari e t al.|2012) . They acknowledge the im- 



portance of S -corrections when doing SNe la photometry with 
multiple instrumental setups: in addition to accurate photomet- 
ric calibration sequences ( |Henden et al.|2012| ), M13 implement 
a light curve merging method that uses phase-dependent zero 
points for each telescope, found via a global x 1 minimization 
per band. For the comparisons with the RqIc photometry, which 
is different from that synthesized using BM12's RI, we used the 
|Landolt|(|1992 ) passbands adjusted to photon count transmission, 



in accord with the standard star catalogs used for calibration by 
|Henden et al.| ( [20T2] l. M13 do not publish individual measure- 
ment errors, but present instead the calibration error for each 
band. This error is smaller than the merged light curve individ- 
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ual night scatter, when multiple observations are available. We 
thus average all same night observations, and use as a conserva- 
tive error the error on the mean of those same observations added 
quadratically to the corresponding band calibration error. The er- 
ror for single observation nights is set to the average of the errors 
of all multiple observation nights. The B-band bias and scatter 
with relation to snifs are the lowest of the three photometric data 
sets considered, attesting to the quality of their calibration pro- 
cedure: M13 is a very good reference for all the B data sets not 
affected by S -correction problems. The biases in the other bands 
are comparable and do not exceed ~ 0.02 mag. The observa- 
tional scatters are slightly larger for VI probably due to outliers 
in one of the 3 merged telescope data sets. The B-band date of 
maximum and peak magnitude found by M13, for which they do 
not quote any error estimate, are very close to our own (within 
0.2 d and 0.02 mag), while the decline rate is perfectly compati- 
ble within our estimated error. 

The synthetic light curves derived from the spectrophotomet- 
ric time series presented in this paper were compared with those 
in the literature. Most of the data sets have the smallest residual 
biases in V, and show an intrinsic scatter of 0.02-0.03 mag for 
all bands. The typical scatter between data sets comes mostly 
from calibration issues, and is of the order of 0.03-0.06 mag, 
similar to the SALT2 fit residuals and with no difference when 
the comparison is made between purely photometric data sets 
or with snifs. This highlights a major benefit of the SNfac- 
tory data, which is at the same time impervious to S -correction 
problems due to its spectophotometric nature, and of photomet- 
ric quality comparable to traditional photometric followup data 
sets, even under the extreme observational conditions pertaining 
toSN2011fe. 



3.2. interstellar absorption in M101 

Estimates of reddening by the host galaxy may be obtained from 
either photometric or spectral data. We consider our photomet- 
ric d^fu^y^i^iiietarijsora 



by |Lira| ( [T996| l and Philli ps et al 



reapplied the procedure used 
( | 1999[ > and derived an intrin- 
sic color law for the "tail" of SN la light curves, and for maxi- 
mum light using pseudocolors (the difference between the mag- 
nitudes of two bands at each band's date of maximum). Their 
color law uses measurements on the CSP photometric system, 
which we synthesize from our spectra using the filters and zero- 
points given by Stritzinger et al. (2011 1. The light curves ob- 



tained are interpolated using SALT2 and two independent fits on 
BVqsp- We measure Am^ficsp) = 1-056 ± 0.035 and obtain 
E(B - V) max = 0.002 + 0.062 mag, by applying Equation 3 of 
|Folatelli et~aT1 ( |2010[ >. We added 0.060 mag in q uadrature to the 
error, to account for the dispersion of the fit by |Folatelli et al.| 
(2010). Applying their Equation 2 to the three spectra whose 
phase is + 30 < t < + 80 d with respect to Vcsp maximum, we 
find a weighted mean and standard deviation of E(B - V) ta ii = 
0.038 ± 0.045 mag. Each individual measurement error was in- 
creased in quadrature by the dispersion of the corresponding fit, 
0.077 mag. Both values are compatible and agree with the trend 
found by Fola telli et al.| ( |2010| l between the difference of both 
measurements and E(B- V) max (upper-left panel of their Fig. 12). 
Using a weighted mean of these estimates, and assuming they are 
independent, we find (E(B - V)host) = 0.026 + 0.036 mag for the 
reddening due to dust on the line of sight to SN 201 lfe in M101. 
This is compatible with the value (0.03 ± 0.06 mag) found by 
Tam mann & Reindl] ( 20 1 1 )l using SN la intrinsic colors derived 
empirically by |Reindl et al.| l |2005) . 



Spectroscopically, we measured simultaneously the equiva- 
lent width of Nai D for the Milky Way (MW) and M101 using 
fixed gaussian profiles with a doublet ratio of 2 : 1 . We obtain 
respectively EW(Na i D) MV 



10+22 m A and EW(Nai D) host 



8 + ii? mA, with a 95% confidence limit at 162 mA. The latter 



-46 



is consistent with the value found by |Nugent et al. (201 lb I us- 
ing HiRES (45 ± 9 mA) and both results are also compatible 
with high resolution measurements by Patat et al. ( 201 l| l, who 
find 38 + 5 mA and 47 ± 2 mA respectively. Using the em- 



pirical relation proposed by Poznanski et al. (20121 to derive 
the dust extinction from the Di + D2 Nai D lines, we obtain 
E(B - V)host = 0.014 ± 0.003 mag, in accord with our photomet- 
rically derived value. Here we find only a statistical uncertainty, 
as the systematic errors are difficult to quantify. 

In conclusion, we confirm there is little evidence for signif- 
icant extinction of SN 201 lfe by dust in its host galaxy. There- 
fore, no reddening corrections are performed on the snifs spec- 
tra, other than for MW extinction. 



3.3. Bolometric light curve and 56 Ni mass 

To construct a bolometric light curve of SN 201 lfe, we combine 
snifs optical spectrophotometry with ultraviolet (UV) and near- 
infrared (NIR) SN la spectral templates. 

The UV template (1600-3400 A) was constructed from 25 
HST STIS spectroscopic observations of SN 201 lfe, encom- 
passing phases - 15 < t < + 24 d. All observations are part 



of the proposal GO- 12298 (Ellis 2009 ) and are publicly available 
from the HST archive. All spectra were taken using the 52x0.2" 
long-slit aperture and one of the UV (G230L or G230LB) or the 
G430L gratings, together covering ~ 1600-5700 A. Same-night 
observations were averaged, and spectra from different gratings 
merged together by averaging over their common wavelength 
range. The ensemble was cut between 1700-3400 A and set 
to decrease linearly to zero flux at 1600 A. Reduction artifacts 
identified as big isolated "spikes" in the spectra were removed. 
A multi-phase SN 201 lfe UV template was then created, by 
linear extrapolation of the time evolution of the flux per wave- 
length bin, for the original 25 phases plus a null spectrum at 
- 20 d, to represent the pre-explosion phase. For each phase of 
SN 201 lfe we intend to reproduce, the template is (achromati- 
cally) flux normalized with respect to the snifs spectrum at their 
common wavelength range (3300-3400 A), in order to account 
for miscalibrations of the HST data. The median of the residuals 
of synthetic photometry performed on the final template, with 
respect to the measurements reported by Brown et al. (2012) 
for Swift/UVOT uvw2, uvm2 and uvwl, are respectively: -0.22, 
-0.21 and -0.20 mag. These values are similar to what was 



found in 3.1.1 once again implying that the UVOT SN 201 lfe 
measurements suffer from large biases. These residuals should 
be taken with caution however, since they are very dependent 
on the blueward flux distribution of the HST spectra, where the 
noise is high and our template is therefore less accurate. Nev- 
ertheless, since the UV flux of a normal SNe la represents only 
a small fraction of the total flux (cf. Fig. the impact of us- 
ing a crude UV template for the bolometric flux will be modest. 
The variance of the template is constructed so that we have a 
0.03 mag error on the full integrated UV flux. 



For the NIR (9700-24000 A), the |Hsiao et al.| ( |2007 ) SN la 
spectral template was used as a starting point. After extrapo- 
lation to each phase, the template is chromatically warped in 
order to match the SN 201 lfe JHK observations of Mathe- 
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Fig. 4. Phase evolution of the ratio of flux in the ultraviolet (1600- 
3400 A) and near-infrared (9700-24000 A) to total bolometric flux for 
SN 201 lfe. The dashed lines represent cubic spline fits. 



son et al. (2012 > from WIYN/WHIRC. A quadratic spline, con- 



structed from the ratio of observed over synthetic photometry for 
each of the effective wavelengths of the three WHIRC filters, is 
used. The variance is constructed from a flat error estimation 
and warped chromatically in a similar way, in order to match the 
Math eson et al.| ( |2~012| ) measurement errors. The median residu- 
als and scatter of the final NIR template with relation to obser- 
vations are at or below the half percent level. 

The bolometric flux is obtained by integrating the full wave- 
length range 1600-24000 A of the UV + snifs + NIR spectra. 
The evolution of the ratio of UV and NIR flux to total flux is 
shown in Fig. [4] The UV flux accounts for a few percent of the 
total flux with a maximum contribution of about 13% attained 
five days before B-band maximum. It then decreases steadily to 
reach about 2% at t > + 20 d. The NIR contribution starts at 
about 15% at very early phases and declines to 5% around 10 
days after B maximum, increasing once again to reach 15-20% 
at the time of the secondary maximum in the NIR (t ~ 30 d), that 
our data did not sample. 

In order to study the time evolution of integrated filter pho- 
tometry relative to the bolometric flux, we simulated the ratio 
to total flux of individual or combined optical passbands, us- 
ing the BM12 throughputs. For the full phase range considered 
here, the individual band with the least scatter on the ratio of 
bolometric flux is /, even if the mean ratio to total flux is small 
(7.3% + 2%). The mean flux ratio on V is higher while keeping 
a relatively small scatter (15.7% ± 3.1%), similar to what was 
seen by | Wang et aL] pool?] ) for SN 2005cf. If we sum fluxes of 
multiple bands, the most promising combinations are those in- 
cluding bands from both extremes of the optical window, namely 
UVI (36.7% ± 1.3%) and URI (38.0% ± 1.5%). The scatter of 
the URI combination decreases below the percent level if one 
limits phases to t < 15 d. The best non-U combination is BI, 
with mean ratio and scatter of 32.1% ± 1.9%. Assuming the rel- 
ative normality of SN 201 lfe, these ratios can thus be used to 
compute estimates of total bolometric flux from near maximum 
optical measurements of a typical SN la. 

The distance-corrected (Shappee & Stanek 201 1 ) bolometric 
luminosity is presented in the last column of Table [2] and the 
bolometric flux is plotted in Fig. [5] along with the integrated flux 
over the full snifs optical window. A f 2 power law was fitted to 
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Fig. 5. Bolometric light curve for SN 201 lfe (filled red circles). The 
open circles represent the integrated flux over the full snifs wavelength 
range and the dashed line is a cubic spline fit for determination of the 
date of maximum. The inset is a zoom of the area marked by a dashed 
box, where a quadratic rise is fit to the first 3 nights of snifs observa- 
tions. The vertical dotted lines mark the time of observations presented 
by |NugenteTaT1P0TTb) . 



our first three observations (inset of Fig. [5]), under the assumption 
that the luminosity dependence is dominated by the evolution 
of the photosphere surface area ( Riess et aL"1|1999[ [Goldhaber 
|etal.|2001[ |Gonzalez-Gaitan et al.|2012[ |Piro|20T2| i until about 
five days after the explosion. We find an explosion date (MJD) 
to = 55796.81 +0.13, where the ~ 3h precision is due to the lack 
of earlier snifs observations. The value is compatible with the 
one found by |Nugent et~aT1 ( |2011b| l (55796.696 ± 0.003) using 
gr-band observations which start two days before the SNfactory 
ones. The goodness of the fit decreases if we include the fourth 



night of snifs observations, supporting the hypothesis that the t 
model is only valid until ~ 5 d after the explosion. By letting the 
exponent of the power law float we find a best-fit exponent of 
2.21 + 0.51 and t = 55796.47 + 0.83. This is a fully constrained 
fit with 3 parameters on 3 points, and earlier observations would 
be needed for better precision. The former (f 2 ) analysis made 
using the Vbmi2 band yields to = 55796.68 + 0.12, in closer 
accordance with the PTF value, the one found by |Brown et al.| 
l|20l2) when using the v band of UVOT (55796.62+0.03), and by 



Vinko et al. 



(2012) (55796.70 + 0.16) using R -band observations. 



The derived data of explosion thus seems to be dependent on the 
wavelengths that are included in the bolometric flux. 

By fitting a cubic spline to the bolometric light curve in 
the phase range - 9 < t < + 15 d, we find the date of max- 
imum bolometric luminosity and derive a bolometric rise time 
T r,boi - tmax,boi ~ k) - 16.58 + 0.14 d, where we use the explo- 
sion date found with the bolometric flux and the f 2 power law, 
and assume an error on the determination of the date of maxi- 
mum of 0.06 d, as with the SALT2 fits. Using the date of Z?-band 
maximum from Sect. 12.41 we obtain t, % b = 17.70 + 0.14 d, 1.1 d 



more than the bolometric rise time and compatible with the SN la 
sample of Contard o et al.| ( [2~000| l, who found both rise-times to 
be within 2 days. One should notice that those authors used 
pseudo-bolometric light curves from integrated UBVRI, most 
closely matched by our snifs integrated flux, from which we find 
a date of maximum 0.5 d later than r r , boi an d hence closer to T r , b- 
Employing similarly integrated snifs data, Sca lzo et al.| ( |2012) l 
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found a difference of ~ 1 d for a sample of overluminous SNe la. 
Our stretch-corrected r' rB — 18.27 + 0.14 d is in agreement with 



the median value found by Gonzalez-Gaitan et al.| ( 2012) for the 
SNLS sample when using SALT2 (18.16 ± 0.44 d). 

The maximum bolometric luminosity of SN 201 lfe obtained 
from the spline fit is L maXiijo / = (1.17 ± 0.04) x 10 43 erg s _1 . 
This can be used along with the light curve rise-time to compute 
the 56 Ni mass synthesized in the explosion. Different authors 
used different rise-time estimates, either from (pseudo) bolo- 
metric or integrated band light curves. As already shown, these 
estimates will differ slightly, thus impacting the computed 56 Ni 
mass. We chose to use our bolometric rise time T r ^oi- Using 
Gonzalez-Gaitan et al.| ( |2012| l, which re- 



Equations 3 & 4 of 



produce Howell et a 
find M 5 6 Ni = (0.44 : 



( |2009) and use y = a = 1.2 ± 0.2, we 
0.08) X {1.21 a) Mo. The a dependency is 
stated explicitly due to its large influence on the final 56 Ni mass 
( |Scalzo et al. 2010[ ), and its uncertainty was propagated in the 
calculations. If one assumes or = 1.0 ± 0.2 (Arnett's rule) as in 
IStritzinger etaL[ ( [2006) , |Wang et al.| p009) and |Hayden et aL[ 
( [2010) , we obtain M 56Ni = 0.53 + 0.11 M G . This value is com- 
patible with the amount of 56 Ni mass that the explosion models 
used by Ropke et aL[(|2012) were set up to produce (~ 0.6 M ). 



3.4. Spectral indicators 

Several spectral indicators defined using maximum light spectral 
features were measured for the SN 201 lfe spectrum closest to B 
maximum light: the 7?si depth ratio and the 7?ca and 7? sis flux 
ratios (|Nugent et al.|1995| |Bongard et al.|2006[>; the ^642/443 flux 
ratio ( [Bailey et al.||2009); the Can H&K (| Walker et al.||2010) 
Si ii ,14131 ([Bronder et al.|2008||Arsenijevic et al.|2008[|Chotard| 



etal.12011) Si H /15972 & ,16355 ( [Hachinger et al.|2006t|Branch1 
et al.f2006||2"009] l and C n A 6580 {cf. Sect.|3 6) equival ent widths. 
The S ii ,15 640, Sin ^6355 ( |Benetti et al.|2004|[2003t|Hachinger| 
|et al.|2006| ) and C n ,16580 feature velocities were also measured 
for all of the spectra for which the features were detectable. For 
most of these spectral indicators (all except "Rsis an d ^642/443X 
a precise estimate of the extrema wavelength and flux positions 
enclosing the feature or defining the velocity is needed in order 
to make a proper measurement. This was done by automatically 
measuring the maximum/minimum flux of the smoothed spec- 
tra in a given wavelength range, and confirming each measure- 
ment visually. The uncertainties on each of these values were 
derived using a Monte Carlo procedure that takes into account 
the impact of the method used to select the feature boundaries, 
the extrema wavelength and flux positions, as well as the small 
contribution of the photon noise. A detailed description of the 
measurement method with the corresponding results on a large 
part of the SNfactory data set will be presented by N. Chotard et 
al. (in preparation). 

The values of these measurements made on the spectrum 
closest to maximum light appear in Table [3] while the phase evo- 
lution of the expansion velocity of S n, Si n, C n and the pseudo- 
equivalent width of C ii ,16580 are plott ed in Fig. [6] T he velocity 
gradient of Si n /16355, as defined by |Benetti et al.| ([2005 ), was 
computed using a linear least squares fit for -2<f<+25d 
overplotted on Fig. [6j and we find u = 59.6 ± 3.2 kms -1 d" 1 . 
The velocity evolution of S n ,15640, which |Benetti et al.|p004) l 
suggests is an effective probe of the true photospheric ve- 
locity, is similar to that observed by those authors for SNe 
1998bu, 1994D and 1990N. It decreases from ~ 12 000 kms" 1 at 
- 15 d to approximately 8 000 km s" 1 at + 10 d, passing through 
~ 9 500 km s" 1 at maximum. 
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Fig. 6. Phase evolution of the expansion velocities of Sn ,15640, 
Sin ,16355, Cn ,16580 (top) and the pseudo-equivalent width of 
Cn ,16580 (bottom). A linear fit for the post-maximum velocity gra- 
dient of Si ii ,16355 is overplotted as a red line. 



Using these spectral indicators and the light curve decline 
rate {cf. Table [3} we place S N 20 life within the SN la clas- 
sification schemes proposed by |BenettTe t al. (2003} and |Branch| 
et al. ( 2006) 1, as seen in Fig. [7] From this comparison, SN 201 lfe 
seems to be a spectroscopically "core normal" SN la whose ex- 



pansion velocity rate of change lies close to the separation be- 
tw een the "low" and "h igh" velocity gradient groups as defined 
by Benetti et al. (2005 v — 70 kms" 1 d" 1 ). It is also a "normal" 
SN la according to the [Wang et"aL] ([2009) definition, based on 
the Sin /16533 expansion velocity at maximum {cf. top panel of 
Fig.§. 

3.5. SYNAPPS fitting 

We analyze the spectral time series of SN 201 lfe with the help 
of the highly parameterized SYNAPPS fitting code (Thomas et al. 
|201 lb) . The objective is to identify ions by their spectroscopic 
signatures taking line-blending into account. The main radiative 
transfer assumptions underlying SYNAPPS are: spherical symme- 
try, a sharply-defined blackbody-emitting pseudo-photosphere, 
pure-resonance line transfer under the Sobolev approximation, 
with line opacity parameterized radially according an exponen- 
tial functional form and per line assuming Boltzmann excitation. 
SYNAPPS combines this parameterized spectrum synthesis cal- 
culation with a parallel non-linear optimization framework to 
reduce the need for tedious interactive adjustment of fit param- 
eters and to assure more systematic sampling of the parameter 
space. The assumptions are simple but sufficient for our imme- 
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Fig. 8. SYNAPPS fits to selected snifs spectrophotometry of SN 201 lfe. Major ion signatures in the first two spectra are shown. The identified 
species are typical of a normal SN la early on and near maximum light. Two weeks after maximum, the spectrum is dominated by numerous 
blends from iron-peak elements (Cm, Fen, Con) blueward of about 5000 A. 



diate purpose; a detailed abundance tomography analysis with 
more sophisticated tools will be the subject of future work. 

Early, at-peak, and two weeks post-peak fits to snifs data ap- 
pear in Fig. [8] Extreme late-time data (after +70 d) are arguably 
more difficult to analyze with SYNAPPS as its physical assump- 
tions are less applicable then. The evolution of the spectral fea- 
tures follows the usual pattern of normal SNe la. Early on, strong 
contributions from O i, Mg n, Si n, S n, and Ca n are detected. 
Cn /16580 and Cn 17234 are detected with high confidence. 
High velocity components for Si n and Ca n features are needed 
to achieve a good fit in the first spectrum. By maximum light 
the high velocity features have weakened and only photospheric- 
velocity opacity components are needed. The C n /16580 feature 
is reproduced by SYNAPPS to some extent at maximum, but the 
fit is not perfect. At about two weeks after maximum, the spec- 
trum is dominated by lines from Fen. We find evidence that 
some Cm and/or Con may be present in the spectrum at this 
point, though [Parrent et aI7| ( |2012) l did not invoke it in their anal- 
ysis. In all other respects we confirm the analysis presented by 
|Parrent et al.1 ( pOT2] ). 

3.6. Carbon signatures 

The spectra of SN 201 lfe contain signatures of unburned car- 
bon from the earliest observations through maximum ( |Cenko| 
|etldT1|2l)TT| |Parrent et aLl[20T2] l. Carbon in SN la ejecta is 
most robustly detected using C n /16580 absorption signatures in 
pre-maximum spectra (|Thomas et al.|2007)|Parrent et al.|2011 



Thomas et al. 201 la| |Folatelli et al.|[20"T2| [Silverman & Filip 
|penko||2012| |Blondin et al.||2012| l. High S/N, mostly nightly 
cadence, and proper treatment of telluric absorptions allow us 



to trace the detailed evolution of relatively weak C n -16580 and 
also -17234 in the snifs time series (Fig. [9]». Both signatures man- 
ifest quite clearly in the earliest spectrum, and gradually fade 
away. C n -16580 is just barely detectable in the +3.7 d spectrum, 
and absent afterward. After -5.3 d, it becomes harder to clearly 
identify the absorption associated with C n 17234. 

In Fig.[6} the blueshift at absorption minimum is plotted. The 
velocity measurements are compared to the ones extracted from 
Si ii 16355 and S n 15640. At all phases, the measured C n ve- 
locity is smaller than that measured from the Si n feature, and 
tracks that of the S n line more closely. This offset could be sen- 
sitive to how strong the Sin line is relative to the other two - 
increasing line opacity in expanding atmospheres (leaving other 
variables fixed) tends to blue-shift absorption minima (e.g. |Jef-| 
|fery & Branch|1990] l. The plot suggests that unburned material 
is present at roughly the same velocities as freshly synthesized 
intermediate mass element ejecta. The carbon appears to extend 
from 15 000 km s _1 (blue absorption edge in the first spectrum) 
down to 8 000 or 9 000 km s" 1 . 

We extract a pseudo-equivalent width for the Cn 16580 ab- 
sorption feature and plot it in Fig.|6]as well. Folatelli et al. (2012 1 
suggest that the pseudo-equivalent width of this feature may in- 
crease and then decrease with time. This behavior was detected 
by |Silverman & Filipperiko] ( |20T2] l in SN 1994D, but the peak 
was not well-sampled. The daily follow-up cadence allows us to 
explore the evolution of this quantity in detail in SN 201 lfe. The 
equivalent width at first decreases for two days, increases back to 
roughly the same value, and then decreases again. This behavior 
was strongly correlated with the velocity evolution in SN 1994D 
(Silverman & Filippenko 2012) but here no such correlation is 
detected. 



Article number, page 13 of 



23 




^ 50 

CM 

| 40 

a 30 

in 

g 20 

LU 

a 10 




□ Broad line 

♦ Cool 

• Core normal 
A Shallow silicon ^ 



♦ 

♦♦ ♦ 



A 



<>1 



□ 



□ 



50 100 150 

pEW (Sill A6355) 



200 



Fig. 7. Top: 7?(Si n) and ii(Si n) versus SN decline rate Am l5 . Compar- 
ison data and sub-classes from Benetti et al. (2005 1, updated with mea- 
surements by Hachinger et al. ( 2006) and Taubenberger et al. (20081. 
Bottom: pEW (Sin A5912) versus pEW ( Sin /i6355). Comparison data 
and sub-classes from Branch et al. ( 2009 ). SN 201 lfe is represented by 
a filled red circle on all panels. 



4. Discussion 

The analysis of the snifs spectrophotometric time series con- 
firms that SN 201 lfe is a "textbook case" SN la. The light 
curve shape and color parameters are close to those of a fidu- 
cial SALT2 SN la (Table[3]summarizes these and other pertinent 
parameters derived for SN 201 lfe in the previous sections). The 
early and near-maximum spectra exhibit typical strong low-to- 
intermediate-mass ion signatures (O i, Mg n, Si n, S n, Ca n) at 
typical SN la ejection velocities. Neither the presence of C n nor 
high-velocity Can absorption (Mazzali et al. 2005) are consid- 
ered particularly unusual. As expected, iron-peak element sig- 
natures dominate at late times as the photosphere recedes deeper 
into the ejecta. Furthermore, there is little evidence for substan- 
tial extinction due to dust in M101 along the line of sight. 

The high cadence of observation, broad wavelength range, 
high S/N, and good calibration make the data set presented here 
not only an asset for studying SN la physics in detail, but also for 
simulations and SN cosmology systematic error analysis. In this 
section, we use our observations of SN 201 lfe to consider a few 
fundamental questions about the conventional analysis of SNe la 
in those areas. This gives us an opportunity to demonstrate some 
generically useful features of the data. 




6000 6500 7000 

Rest Frame Wavelength (A) 

Fig. 9. Evolution of Cn /16580 and ,17234 features in the spectral 
time series of SN 201 lfe. Two axes depicting blueshift with respect 
to 6580 A and 7234 A are overlaid, with velocities given in units of 
10 3 km s -1 . Cn ,16580 begins as a notch which gradually recedes in 
velocity and weakens progressively until it is difficult or impossible to 
reliably identify at + 6.7 d. C n /17234 follows a similar evolution, but 
disappears by - 1.3 d. 



4.1. Comparisons with other SNe la 

Figure [10] compares selected snifs observations of SN 201 lfe 
with some of SNF200805 14-0 02 ( |Aldering et U] [2008) ) and 
SN 2009ig (Kleiser et al. 2009]). Here we normalize the mean 
flux between 3500 and 7500 A in each spectrum involved, and 
overlay matched phases. The relative flux as a function of phase 
is not preserved, but this allows us to compare spectral features. 
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Fig. 10. Comparison of selected epochs of SN 201 lfe with SNF200805 14-002 (left, panel a) and SN 2009ig (right, panel b). All observations 
were obtained with snifs. Insets in each panel show the region around the Cn /16580 absorption notch (5800-6800 A rest-frame) in the earliest 
pair of spectra compared. 



These two SNe la are selected for comparison because they have 
been identified previously as basically "normal" SNe la ( [Thomas | 
et al.|201 la[ [Foley et al.|20T2a) . 



SNF200805 14-002 and SN 201 lfe are a good match at all 
phases compared, but they are not identical twins (Fig. [TOfr). 
The agreement is best redward of 4000 A, in terms of what fea- 
tures are present and their morphologies. For example, both ob- 
jects have strong C n /16580 absorption notches early on (see in- 
set, near 6300 A). The double notch S n features near 4700 A 
and other small-scale features are also matched well. In con- 
trast, however, the - 10.2 d spectrum of SN 201 lfe exhibits 
a very slight high- velocity Can infrared triplet absorption that 
is weaker or absent from the corresponding SNF200805 14-002 
spectrum. This is mirrored in the near-UV where the Ca n H&K 
absorption feature extends to higher velocity in SN 201 lfe than 
in SNF200805 14-002. In fact, at this phase the spectra blue- 
ward of 3800 A are markedly dissimilar. Five days later, the 
morphologies of the near-UV spectra appear to have converged 
together. About three weeks past maximum, some overall color 
difference is apparent between the two SNe, but again the spec- 
tral features are quite similar. 

There is substantial evidence in the literature for greater di- 
versity in the near-UV properties of SNe la than in the optical 



( Ellis et al 120081 |Brown et al.(20T0l |Milne et al.|2010| |Cooke 



et al.|201 l[|Wang et al.|2012[|Foley et al.|2012b| >. Thenear-UV 
behavior depicted in Fig.|10p supports the idea that this extends 
to the time domain as well. Increased diversity in the near-UV at 
early times relative to later phases seems plausible from a phys- 
ical standpoint. This region of the spectrum is particularly sen- 
sitive to the composition, density, and temperature of the outer 
layers of the ejecta, and at early times the spectrum forms in 
these layers. 

SN 2009ig (Fig. [TO]? ) provides an excellent foil for 
SN 201 lfe. |Foley et al] ( |2012a| > noted the relatively high 
blueshifts of absorption features in its spectra, and here we can 
see the contrast with the more typical blueshifts exhibited by 
SN 201 lfe. In the earliest spectrum we clearly see the system- 
atically higher velocity absorptions of Ca n H&K and infrared 
triplet and Sin A6355 in SN 2009ig. The blue edges of these 
features weaken with time, but overall the absorption lines re- 
main faster in SN 2009ig than in SN 201 lfe. An inset compares 



the region of the spectrum around the C n /16580 notch in the 
first spectrum of either SN. A robust notch is seen in SN 201 lfe, 
but any Cn absorption in SN 2009ig at the same phase is much 
weaker ( |Parrent et al.|201 l[|FoTey et al.|2012al >. 

The above meta-comparison demonstrates a key benefit of 
the dense temporal sampling of our SN 201 lfe spectrophotome- 
try, that it is a trivial matter to find SN 201 lfe spectra to compare 
to spectra of other SNe at practically exactly the same phase. 
This should be a highly useful feature of the data set for other 
researchers who wish to contrast their observations with a clas- 
sic, normal SN la in a systematic way. 

4.2. Comparisons with spectral surface templates 

Here we examine how well the color and detailed spectral evo- 
lution of SN 201 lfe are described by spectral surface (A, t) 
templates used for conventional light curve fitting and K- and 
S-corrections. Such templates should do well at describing 
"archetypical" SNe la. Our default expectation is that the tem- 
plate spectra should have high fidelity to the "ground truth" time 
series, that of SN 201 lfe. The comparison allows us to gauge 
the level of that fidelity, and identify transient or systematic defi- 
ciencies in the templates themselves. We begin by examining the 
agreement between true and template-synthesized color curves 
and then take a closer look at the agreement at the fine spectral 
level. 

B - V, V -R, and V - I color curves of SN 201 lfe, syn- 
thesized from snifs spectrophotometry (Sect. 2.3 i up to 40 days 
past explosion, using BM12 passbands and zero points, appear in 
Fig.fTTj Color cu rves synthesized from standard templates (Nu- 
gent et al.|2002| |Hsiao et al.|2007[ ), corrected to the stretch and 
observed colors at peak luminosity, are shown for comparison. 

The observed B - V color curve closely follows the Hsiao 
template starting at about - 10 d, while the match with the Nu- 
gent template is only good after about - 5 d and diverges slightly 
starting at +5 d. For very early epochs (before - 10 d) neither 
template seems to be a good match. 

In the case of the V - R color evolution, SN 201 lfe matches 
reasonably well with the Hsiao template after maximum light, 
while earlier phases differ significantly. Interestingly, the Nugent 
template shows a large deviation from the Hsiao template at the 
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Table 3. Relevant parameters for SN 201 lfe derived from this work 



Parameter 


Value 


Photometry 


Xi 


-0.206 + 0.071 


Color 


-0.066 ± 0.021 


Stretch a 


0.969 ± 0.010 


Am i 5 B h 


1.103 ± 0.035 


^max *J 


55813.13 + 0.06 


^max $ 


55814.51 +0.06 


Wx v 


55816.25 + 0.06 


^max ^ 


55816.06 + 0.06 


^max ' 


55812.55 + 0.06 


^ max / *-* maxg 


9.49 + 0.02 / 9.52 + 0.02 mag 


^max 


9.94 + 0.01 mag 


V / V 

v max / " maxg 


o os + 02 / 9 99 + 01 mas 


^max / ^maxg 


10.02 + 0.04 / 10.04 + 0.03 mag 


Anax / Anaxs 


10.30 + 0.11 / 10.32 + 0.09 mag 


M max 17 


-19.55 ±0.19 mag 


M max fi 


-19.10 + 0.19 mag 


M max y 


-19.06 ± 0.19 mag 


JVl max K 


1 o no j. n iq 
— ly.yjz + \j.iy mag 


M / 

max 


-18.74 + 0.22 mag 




0?6 + 036 


e(b - y)host spetlI0 


0.014 + 0.003 


fo bolometric d 


55796.81 +0.13 


to V d 


55796.68 + 0.12 


r r bolometric 


16.58 ±0.14 d 


T r B 


17.70 ± 0.14 d 


j max 

Hoi 


(1.17 ± 0.04) x 10 43 ergs s" 1 


M56 Ni 


(0.44 ± 0.08) x(1.2/a)M 


Spectroscopy 




0.26 




1.27 


7?SiS 


1.32 


^642/443 


0.73 


pEW (CanH&K) 


97.62 ± 11.35 A 


pEW(Sin/l4131) 


15.33 ±0.07 A 


pEW (Sin ,15972) 


17.21 ±2.19 A 


pEW (Sin ,16355) 


96.67 ± 0.53 A 


pEW(Cn ,16580) 


0.56 ±0.17 A 


v (Si ii ,16355) 


59.6±3.2kms-'d- 1 



Notes. All dates are MJD. Magnitudes are on the "standard" SALT2 
magnitude system, allowing direct comparison with the literature. Ab- 
solute magnitudes are computed using the distance modulus to M101 
given by Shappee & Stanek (201 1 1. All spectral flux ratios and pseudo- 
equivalent widths are derived at maximum light. The photon noise de- 
rived uncertainties for the flux r atios are below th e percent level. 
Derived using Equation 6 of Guy et al7| < |2010^ . 

(b) Derived using the fitted SALT2 light curve model. 

(c) Observed magnitude at the time of S-band maximum. 

(d) Derived assuming a t 2 power law. 



same phases that the Hsiao template differs from the SN 201 lfe 
observations; notably the "bumps" at - 10 d and (to a lesser ex- 
tent) + 5 d. While this difference is not as pronounced as the 
one between both templates, one could argue that in the small 




sample of SNe la used to construct the Nugent et al. 



plate, there are spectra whose time evolution is c 



(2002) tem- 



oser to that 



-10 -5 5 10 15 20 
Phase Relative to B Maximum (days) 

Fig. 11. Synthetic color curves of SN 201 lfe from snifs spectropho- 
tometry. Both spectral templates from Nuge nt et al.| |2002|; |Hsiao| 
|et al.|(2"0"07| l are corrected to the observed stretch and peak colors. The 
SALT2 curves are derived from the fitted light curves on the different 
bands. The error bars on the SN 201 lfe values are smaller than the plot 
symbols. 



of SN 201 lfe than the average SN la spectral time series rep- 
resented by the |Hsiao et al.| ( |2007] > template. This appears to 
be confirmed by the fact that the two SNe la that SN 201 lfe is 
most similar to spectroscopically (SN 1992A and SN 1994D, ac- 
cording to Nu gent et al.|201 lb I are those contributing the most 
spectra to the Nugent template, along with SN 1989B. 

The evolution of V - 1 is similar to that of V - R. After max- 
imum light, the data agree rather well with the Hsiao template, 
with the Nugent template being a better match from + 15 d on- 
wards. At earlier phases the observations display a significant 
systematic departure from any of the templates, being closer to 
the Nugent template at about - 10 d. From both this color curve 
and the previous, it is apparent that the second snifs spectrum 
(MJD 55800) has a problem at redder wavelengths, being fainter 
than we would expect it to be. This systematic error (~ -0.03 
and ~ -0.04 mag for RI respectively) is close to the statistical 
uncertainty for the flux calibration of that night and seems to be 
due to extraction problems, not with the SN itself but rather with 
the sole standard star used for calibration. 

The "fitted" SALT2 color curves and error bands are also 
plotted in Fig.[TT] For these, we can see that the fitted colors are 
accurate throughout all of the phases after - 5 d, but start to break 
away from the observations at phases $ - 10 d for B—V and be- 
fore - 5 d for both color curves using the redder parts of the 
spectra. For the latter, the fits are systematically too red, though 
the observations are within the quoted model errors. This is not 
unexpected since the S ALT2 model is trained on both nearby and 



Article number, page 16 of 



23 



R. Pereira & the SNfactory: Spectrophotometric time series of SN 201 lfe 



high-z SNe la observations, with the former often lacking obser- 
vations at very early phases, and the latter being affected by low 
signal-to-noise for the redder bands due to quantum efficiency 
drops and atmospheric extinction effects. 



Next we consider agreement at the fine spectral level. Fig. 12 



is a comparison of six representative snifs observations of 
SN 201 lfe with two SN la spectral surface templates available 
in the literature, that from |Hsiao et al.| ( |2007| l and SALT2. The 
rapidly changing spectral features during the rise in luminosity to 
peak are the most interesting to consider, so four spectra leading 
up to maximum are shown. Observations after peak brightness 
are represented by the spectra at + 1 1.7 and + 23.7 d. 

The Hsiao template spectra were warped as in |Foley et al.| 
( |2012a[ l, in order to emulate their usage in K- or 5-correction 
computations. For each phase the template spectrum is normal- 
ized by a cubic spline fitted to the five ratios of synthetic UB- 
VRIsm fluxes between the template and our spectrum for that 
same phase. In this way the template is forced to match the ob- 
served broad-band colors of SN 201 lfe. 

The SALT2 spectra, on the other hand, are created directly 
from the template using the light curve fit parameters (xq, X\ 
and c) found in Sect. 2.4 without any additional warping over 
wavelength or time. Thus, the projected spectrum from SALT2 
is what it predicts the underlying SN spectrum looks like with- 
out modification. No scaling in flux is added to improve the 
agreement, even cosmetically. The SALT2 model lcr error is de- 
picted as the gray shaded region in each panel of Fig. [12] The 
lower plot in each panel shows the percent residual differences 
between each template and the observations. 

Though there are both transient and persistent artifacts in the 
residuals, a general trend is also evident. The agreement be- 
tween the templates and SN 201 lfe improves as maximum light 
is approached. Near maximum, where the population of SNe la 
available for constructing either template is very well sampled, 
the discrepancies are relatively smaller. After this point, glitches 
in the residuals become larger again. We note that the residu- 
als are decidedly not the result of a bias induced by comparing 
the spectra to templates scaled by synthetic photometry, the high 
S/N of the observations limits this effect to less than 1%. 

We start by examining the agreement between SN 201 lfe 
and both templates at phases t < -10 d. The disagreement be- 
tween the projected SALT2 spectrum in Fig.[T2]and SN 201 lfe 
at -15.2 d is at first striking. The SALT2 error model strongly 
deweights this observation (and others at the earliest phases, 
see Sect. |2.4) , so the mismatch is not ultimately catastrophic. 
To an extent, this suggests that these early epochs are under- 
utilized in modeling SNe la for cosmological applications to- 
day. By -12.2 d, the template shows signs of converging toward 
SN 201 lfe but both the estimated errors and residuals remain 
large. The Hsiao template is not a particularly good match to 
SN 201 lfe at these phases either. In particular, it seems to sys- 
tematically underestimate the depths and widths of the stronger 
absorption features. The template is based on a comparatively 
small number of spectra at early phases: 13 spectra with phases 
between -15 and -10 d, compared to 50 spectra or more in each 
5 -day bin from -10 to +15 d. The template is obviously much 
more susceptible to biases in the input sample before -10 d. 

SALT2 systematically under-estimates the flux at wave- 
lengths redward of about 7500 A, at all phases. Setting aside 
the case of -15.2 d, the SALT2 -predicted flux at these wave- 
lengths always appears to be smaller than the true value, and 
the shape of the broad Ca n infrared triplet feature is poorly re- 
constructed. The warped Hsiao template seems to have higher 
fidelity to SN 201 lfe at these wavelengths than SALT2 in all but 



the first two spectra shown. While there are similarly sized off- 
sets at the bluest wavelengths, they are not either systematically 
low or high as at the red end. 

Interestingly, at all phases where the bellwether Si n /16355 
absorption is unblended and easily discernable in the observed 
spectrum, it corresponds to a prominent, isolated, coherent fea- 
ture of the residuals. The same applies to the Can H&K feature. 
Discrepancies in feature strength, width, and blueshift are all to 
blame. Even at maximum light, the residuals relative to SALT2 
are larger than the model error in certain places (by factors of 
several at some wavelengths) while the integral of those residu- 
als in the photometric band considered remains small. This em- 
phasizes the difficulty in extracting accurate spectral templates 
by combining spectra with photometry, and the usefulness of 
spectrophotometric SN follow-up for this purpose. In the last 
spectrum, we see that most of the features are reconstructed at 
the right wavelength, but the emission profiles in either template 
are too low or too high. 

Considering the SALT2 residuals, it at first seems plausible 
that discrepancies such as those noted here are symptoms of in- 
completeness in the sample used to create the template. Merely 
expanding the sample would then ameliorate any resulting er- 
rors. This suggestion is certainly applicable to the earliest ob- 
servations, but seems less reasonable closer to maximum light. 
The fact that the SALT2 template spectrum at -0.3 d very closely 
resembles that of a normal SN la suggests that quantity and qual- 
ity of the underlying sample are not the only issues. Rather, it 
seems plausible that additional information in the form of one 
or more model parameters (beyond light curve shape and color) 
are needed to null out coherent residuals associated with major 
spectral features, in particular those originating in the Ca n H&K 
and Si n /I6355 features (see also Cho tard et al.|201 1) . It seems 
quite plausible that in these features the residuals arise from ve- 
locity mismatch, and that if spectra were used to constrain the 
template fit, the contribution to the photometric correction error 
could be reduced. Of course, spectrophotometric measurements 
eliminate the problem entirely. 

It is important to acknowledge that comparing templates de- 
signed to stand in for archetypical SNe la to SN 201 lfe may 
suggest an optimistic outlook in the context of light curve fit 
systematics. On the other hand, it provides a near best-case or at 
least typical scenario. To contrast, we refer the reader to a similar 
comparison of SN 2009ig and the Hsiao et al. ( 2007 1 template by 
|Foley et aT] ( |2012a] ). In that particular case, the discrepancies in 
spectral evolution presented were portrayed as severe. We note 
that we performed the same projected SALT2 residuals analysis 
as depicted in Fig. [12] for SN 2009ig, and find the same kinds 
of residual structures with often larger but at times smaller mag- 
nitude than is seen in SN 201 lfe. A more systematic analysis 
of spectral surface templates is underway and will be presented 
elsewhere (Saunders et al., in preparation). 



4.3. Detecting C n in SN la spectra 

Since SNe la are likely the thermonuclear incineration of carbon- 
oxygen white dwarfs, estimates of the total mass and spatial dis- 
tribution of unprocessed carbon in SN ejecta provide a way to 
constrain the explosion mechanism. Pure turbulent deflagration 
models predict that a large amount of unprocessed carbon may 
remain after explosion (e.g., Gamezo et al. 2003 ; Rop ke & Hille- 1 
brandt 2005 ; Ropke et al. 2007 1. Models where this deflagration 



is followed by a detonation phase consume much of the remain- 
ing carbon (e.g., |Hoflich et al.l|2002] |Kasen et al.||2009] >. Oxy- 
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Fig. 12. Comparison of six different phases of SN 201 lfe (black line) with the Hsiao et al. 1 2007 ) (black dashed lines) and SALT2 (red line 
with grey error model) SNe la spectral surface templates. The Hsiao et al. (2007J spectra are warped to match the observed broad-band colors of 
SN 201 lfe. SALT2 spectra are created directly from the template using the light curve fit parameters. The lower panels plot the percent residuals 
of each template with relation to the observed SN 201 lfe spectra. 
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gen signatures in the photospheric phase are less useful, because 
oxygen is both a fuel and product in the explosion. 

The high quality of the snifs time series and the clear pres- 
ence of C ii 16580 provide us with the opportunity to explore 
some questions about the efficiency of carbon detection in SN la 
spectra. SN 201 lfe gives us a "ground truth" Cn 16580 de- 
tection at high S/N needed to simulate the impact of noise and 
spectrograph resolution on this principal carbon indicator. The 
feature manifests as a small notch or small-amplitude depression 
in flux, so both noise and resolution are important factors in its 
(positive) detection. |Parrent et al.| ( |20lT] l briefly explore the ef- 
fect of noise, but only on a single illustrative example spectrum. 
Here we perform simulations at several phases, especially later 
than - 15 d which is more typical for initial SN la follow-up. 
Specifically, we are interested in simulating quick analyses done 
by observers attempting to assess whether C n -16580 is present 
in freshly obtained spectra, as is commonly described in astro- 
nomical circulars. 

Figure [13]demonstrates the effect of lowering the S/N in five 
pre-maximum spectra of SN 201 lfe where C n -16580 is clearly 
evident in the snifs data. The real observations of the C n -16580 
region, adopted as effectively perfect, appear in the top row of 
panels in the figure. The same spectrum is shown on each sub- 
sequent row with gradually decreasing S/N (given as S/N per 
rest-frame A in the neighborhood of the carbon notch). At a S/N 
of 10 per rest- frame A (bottom row), the C n notch is completely 
obscured at the native snifs extraction binning (2.38 A). Over- 
all, the figure suggests that relatively low S/N may be tolerable 
at early times for making a correct carbon detection, because the 
feature's signal is strongest then. Positive identification becomes 
more difficult as the feature fades away. Indeed, it is hard to 
make a case that the C n 16580 notch is clearly present at -1 .3 d 
at any simulated S/N. 

Binning the data can recover the feature, but at the cost of 
effectively smoothing the data, as shown in Fig. [14] Simulations 
generated as in Fig. 13 are binned by a factor of eight. This is 



particularly detrimental at the earliest phases and, somewhat sur- 
prisingly, at the lowest S/N ratios. The small wavelength scale 
of the feature means that it is susceptible to contrast losses if 
the data are simply binned by too large a factor. We suggest 
that binning may help improve confidence in more marginal C n 
detections at low S/N, but not always. Also, systematic spec- 
tral feature modeling (e.g., Thom as"et al.|201 la| ) may eventually 
provide a quantitative approach. 

Finally, Fig.[l5]depicts simulated data as before, but with the 
resolution degraded to R ~ 100. The top row again shows the 
snifs spectra without any noise added, but with a Gaussian filter 
applied. Only a single noisy realization is plotted here, since it 
is clear from the top row of the figure that at this resolution the 
C ii 16580 feature is obliterated nearly completely at all phases. 
Unless careful modeling is used, it is doubtful that a positive 
carbon detection would be made at any phase. A characteristic 
"break" at about 6250 A at the earliest phase in the simulated 
spectrum may be a useful signature, but could be indistinguish- 
able from normal variation in the morphology of Si ii 16355 fea- 
tures. 

The above reinforces the somewhat obvious fact that con- 
structing an accurate descriptive census of carbon in SN la ejecta 
requires both high S/N and moderately high resolution. Earlier 
is better than later, but blue-shifted C n 16580 notches "blend" 
more readily with the Si n 16355 absorption at lower S/N. Lower 
resolution robotic integral field spectrographs, such as the SED 
Machine (Ben- Ami et al.|20"l2] l, are an excellent way to provide 



prompt broad classifications and phase estimates for SNe, but 
are unlikely to directly contribute to a carbon census in SNe la. 
This will fall to higher resolution spectrographs with larger aper- 
tures that can trigger on early SN la alerts from first-response 
screening that such instruments as the SED Machine provide. 
Real-time forward modeling of features that includes instrumen- 
tal effects may be critical to make the most of such data. 

Given the difficulties inherent in translating early discovery 
into early spectroscopic follow-up of SNe la, it seems quite pos- 
sible that photospheric-velocity carbon signatures at the earliest 
phases may be anything from common to ubiquitous. Consid- 
ering also that C n line strengths vary from object to object and 
as a function of phase (e.g., Parrent et al.||201"T| |Thomas et al.| 
201 la I, the clearest route to a proper census is also the most 



challenging: orchestrating discovery and follow-up on the same 
or next night for day-old SNe la. This remains far from rou- 
tine for wide-field surveys (such as PTF). High-cadence surveys 
targeting nearby galaxies (e.g., the Lick Observatory Supernova 
Search, Filippe nko et al.|2001) > may slowly build up a sample of 
bright objects. For now, it seems wise to consider the current 
rate estimates from the literatur e (|Parrent et al.|2011 |Th omas 



et al.|201 la||FolateUi et al.|2012| [Silverman & Filippenko|2012 



Blondin et al. 2012) to be lower limits. 

The qualitative study presented here highlights the useful- 
ness of high S/N, flux calibrated time series to study detection 
systematics of narrow and weak spectral features. Our data are a 
natural input for systematics study simulations where the accu- 
rate shape of the SN pseudo-continuum is important. Using any 
of the current templates for the same purpose will yield biassed 
results, since they are by nature capturing an average behavior, 
tuned to be useful in a different context. 

5. Conclusion 

We have presented a new atlas of spectrophotometry of 
SN 201 lfe, a photometrically and spectroscopically normal 
SN la. SN 201 lfe is not the typical SNfactory target, and hence 
its calibration route is not the standard one used for the bulk of 
the SNfactory data set. Difficulties encountered were the result 
of short exposure times and very high airmass. Still, the attained 
level of calibration, high S/N, and observing cadence make this 
data set extremely useful for a variety of science applications. 

The SALT2 fit performed on synthesized light curves shows 
SN 201 lfe to have attained a B-band maximum of 9.94 + 
0.01 mag on MJD 55814.51 + 0.06, with "standard" lightcurve 
shape parameters (x x = -0.206 + 0.071, c = -0.066 + 0.021). 
The residual scatter when compared with published SN 201 lfe 
photometry is comparable to the errors estimated for the flux cal- 
ibration, and of the same order as the S ALT2 fit residuals and the 
scatter between independent photometric followup campaigns. 
Reddening due to dust in the host galaxy is found to be very 
moderate and in accordance with independent spectroscopic de- 
terminations. An ultraviolet + optical + near-infrared template 
was built from snifs data and public UV spectroscopy and NIR 
photometry, to construct a bolometric light curve. From it we 
derive a date of explosion fo = 55796.62 + 0.13 in the assump- 
tion of a i 1 luminosity evolution for the very early phases, with 
a rise-time r r = 16.58 + 0.14 d and an inferred 56 Ni mass of 
(0.44 + 0.08) x (1.2/of)M Q . 

The analysis of spectral indicators shows SN 201 lfe to be 
a spectroscopically core-normal SN la, on the lower side of 
the separation between the HVG and LVG groups as defined 
by the rate of change of the expansion velocity of Si ii 16355 
(59.6 ± 3.2 kms~' d _1 ). The evolution of spectral features is 
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Fig. 13. Effect of noise at various phases on the C n .46580 absorption notch. The top row depicts the neighborhood of the feature as observed 
by snifs. Each successive row simulates progressively noisier observations, with the S/N per A labeled at right. The red vertical line depicts the 
center of the C n notch to guide the eye. Early observations are no guarantee of a C n detection even if it is there at lower S/N but generally the 
feature becomes harder to discern with time as the notch amplitude decreases. Wavelength and flux scales depicted are the same in each panel. 
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Fig. 14. Same as Fig. 13 but rebinned by a factor of eight (bin size ~ 19 A). Rebinning noisy data may improve confidence in C n /16580 detection 
in some cases. 
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Fig. 15. Same as Fig.[T3]but with the spectral resolution degraded to R ~ 100. The top row depicts the smoothed input to the simulation, which 
nearly completely destroys the Cn signature at all epochs. At the earliest phase a "break" in the shape of the Sin /16355 feature can be seen near 
6250 A. The break signature may be a way to detect C n at lower spectral resolutions and binning than our observations here, but only at the earliest 
phases. 
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typical of a normal SN la, with low to intermediate-mass el- 
ements dominating the early spectra, and iron-peak signatures 
strengthening after maximum light. High-velocity Ca n and Si n 
are needed to explain the early spectra, but these weaken by 
maximum light. C n /I6580 and C n A7234 are detected in pre- 
maximum spectra, and our high S/N observations and daily ca- 
dence allow us to follow the evolution of these features very 
closely as they fade: apparently unburned carbon extends as 
deep as 8 000 km s in some normal SNe la. 

The snifs time series possesses a number of features that 
should make it highly useful for interpreting other SNe la, inves- 
tigating systematic errors in traditional light curve analysis, and 
making forecasts for new instruments' sensitivities to detecting 
physically important features. 
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(placeholder for animation) 

Fig. .1. Animation of the spectral evolution of SN 201 lfe for phases - 15 < t < + 45 d. The snifs spectrophotometric data was interpolated 
linearly per wavelength bin between each observation date. The inset shows the light curves synthesized using the UBVRIsNf filter set, represented 
as colored shaded regions. 



